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Abstract 

The fractal dimension Df of sites resisting Landau or maximal 
Abelian(MA) gauge fixing in lattice SU(3) gluodynamics is defined 
and computed. In Landau gauge such sites clump into _Dj ~ 1 clus- 
ters in the confining phase. In the finite temperature phase their 
dimensionality drops to _Dj < 1, that is, clustering seems to dissipate. 
In contrast, MA gauge resistant sites fail to exhibit a notable tendency 
to cluster at any temperature. 



Gauge fixing is an essential step for an increasing number of lattice QCD 
applications ranging from the determination of quark wavefunctions to the 
study of the role of magnetic monopoles in quark confinement [1]. As prac- 
titioners know, sites resistant to gauge fixing are routinely encountered in 
Monte Carlo QCD gauge configurations. In this Note, we describe the space- 
time distribution of these sites and show, by measuring their fractal dimen- 
sionality, that resistant sites are not always randomly distributed. Rather, 
in Landau gauge they cluster into string-like formations and, further, the 
fractal dimensionality of these formations is highly temperature sensitive. In 
contrast, resistant sites in another popular gauge, maximal Abelian gauge, 
does not clearly exhibit clustering at any temperature. 

Landau gauge-fixing in lattice gauge theories is achieved by maximizing 
the function 

^ ^ F^{x), F^{x) ^ ReTr(K[/.,,yJ+.) (1) 

x(^A 11=1 

with respect to gauge transformations Vx- Ux,fi refers to the link fields. Since 
F^{x) at any one site can be maximized exactly by suitably choosing V^, 
in the simplest algorithms is maximized iteratively by sweeping through 
the lattice (perhaps in a checkerboard scheme) and maximizing F^{x) at 
each site [2]. However, since nearest neighbor gauge-transformations at x it /i 
ruin the maximization of F^{x), convergence to a global or even a metastable 
maximum of JF^ is not guaranteed. Nonetheless, in practice this "relaxation" 
algorithm does increase JF^ and the gauge configurations do converge to 
something which at most sites looks like d^A^ = 0. 

As previously noted by several authors [2, 3, 4], the relaxation procedure 
is a limiting case of a Monte Carlo simulation of a quenched Higgs model. 
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In this formulation, gauge transformations comprise an adjoint SU{N) 
Higgs field coupled to the gauge fields by the action Sh = —k.J-'^. k is a real 
parameter. The Higgs field is quenched because it also interacts (nonlocally) 
with the gauge fields via the Faddeev-Popov determinant whose inverse is 

A^p(«;{t/.,4) = f[dV] exp{-SH}. (2) 

Landau gauge corresponds to the light Higgs or large k, limit; taking k, be- 
yond some critically large value Kc forces App to be dominated by Higgs 
configurations which maximize JF^. 

At any value of k, A]^p can also be viewed as the partition function of an 
SU{N) spin model with locally variable, pseudo-random couplings {nUx,^}- 
(Only when fiqcD = are the couplings completely random.) From this 
viewpoint [5], relaxation corresponds to doing n > Hc Monte Carlo sweeps 
for the semiclassical ground state spin configurations. Therefore, analogous 
to defects produced in the instantaneous freezing of liquid iron, defects are 
expected to be produced during relaxation. 

Indeed, a histogram of the F^{x) distribution after Landau gauge fixing 
reveals a small but long tail of small-F^(a;) "resistant" sites. As we will 
show, resistant sites are not homogeneously distributed in space but tend to 
cluster into spacetime regions, "defects." It is important to get a quantitative 
handle on defects since they surely affect both action Sh and Faddeev-Popov 
determinant App and, hence, may potentially infiuence lattice results which 
depend on gauge fixing. 

We must emphasize that, as previously illustrated in Ref [4], some defects 
are essential features of some lattice gauge theories and not just artifacts of 
relaxation which may be potentially suppressed by, for example, an anneal- 
ing procedure. In compact QFD essential defects occur because magnetic 
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monopoles — which are gauge invariant and physical in this model — populate 
the gauge configurations. All magnetic monopoles are connected to their 
antimonopole partners by Dirac strings, which are thusly unavoidable (al- 
though gauge variant). Since the Landau gauge vector potential around 
a Dirac string in the z direction is (in continuum cylindrical coordinates) 
^string _ (j)lf^^ the Landau gauge value of F^{x) = X]^=i cos A^(a;) must be 
disrupted near a Dirac string or, equivalently, along a path connecting each 
monopole-antimonopole pair. As depicted in Figure 1 of Ref. [4], resistant 
sites cluster around Dirac strings. In D = 3 + 1 dimensional spacetime, Dirac 
strings sweep out worldsheets which would cause fractal dimension ~ 2 
Landau gauge defects. Such defects are indicative of the disorder caused 
by Dirac strings which, for example, gives the photon propagator a nonzero 
mass pole in the confined phase. is gauge dependent because the lo- 
cation and density of Dirac strings vary with gauge. In the deconfined phase, 
monopoles and Landau gauge Dirac strings become dilute. Correspondingly, 
resistance to gauge fixing dissipates and in Landau gauge vanishes. 

In this Note we show that the fractal dimension of Landau gauge defects 
in the confining phase of SU (3) gluodynamics is also nontrivial. Furthermore, 
in the finite temperature phase we find that these defects tend to dissipate 
and Df is dramatically smaller. 

For comparison, we also look for defects in maximal Abelian(MA) gauge [6], 
putatively used to identify nonAbelian magnetic monopoles in 't Hooft's 
Abelian projection scheme [7]. MA gauge fixing is achieved by maximizing 
-^^ = E.GA F^ix) where 

a /i=l 

and is the sum over the Cartan generators of the gauge group. Our MA 
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gauge defects are defined using F^{x) in complete analogy to Landau gauge 
defects. As described below, resistant sites do not exhibit pronounce clus- 
tering in MA gauge at any temperature. Since both MA and Landau gauge 
are numerically achieved using the same checkerboard relaxation sweeping 
scheme, the absence of defects in MA gauge tends to rule out algorithmic 
artifacts as the cause of clustering in Landau gauge. 

A compelling difference between MA and Landau gauge is that, unlike 
Landau, MA gauge leaves a residual [t^(f)]^~^ local gauge invariance: the 
set of all transformations which commute with the Cartan generators 
in (3). This means J-'^ is insensitive to structures which obstruct gauge 
fixing of the subgroup whereas JF^ is not. Could Landau gauge 

defects be manifestations of such obstructions? Naively, one might 

anticipate the following. Consider the world lines of the aforementioned 
Abelian projection monopoles. Such monopoles, typically identified on the 
lattice by fixing to MA gauge, are invariant and carry 

magnetic charge. Every monopole is connected to its antimonopole partner 
by a variant Dirac string. As in compact QED, these Dirac strings 

would sweep out world sheets with a fractal dimensionality of ~ 2. Hence, 
they might induce ~ 2 defects in Landau gauge fixing which do not appear 
in MA gauge fixing. 

Unfortunately, our simulations indicate that the actual situation is not so 
transparent. As described below. Landau gauge defects have an apparent Dj 
which is closer to f than 2. This discrepancy may rule out Dirac worldsheets 
as the culprits. The only known ~ f structures in QCD are the magnetic 
monopole currents themselves. (Fractal dimensionality of some other QCD 
features is discussed in [8].) Yet, being invariant, there is no 

obvious reason why they should make defects in Landau but not MA gauge. 
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Additionally, as monopole identification varies with the gauge of the Abelian 
projection it is hard to write down a concrete relationship between the MA 
gauge monopole currents and the Landau gauge defects. As described below, 
we have performed numerical spacetime correlation studies which indicate 
that Landau gauge defects are not closely correlated with the locations of 
the MA gauge monopole currents. This negative result does not necessarily 
preclude the possibility that defects and monopoles are manifestations of the 
same underlying QCD effect. It could simply be that these manifestations 
shift locations between Landau and MA gauge. 

The numerical results to be described are based on the following fiqcD = 
6.0 pure gauge 24"^ X Nt SU{3) lattices: eight Nf = 40 configurations, nine 
Nf = 8 configurations, five Nt = 6 configurations, and eight A^t = 4 config- 
urations. Nf is the extent of the Euclidean time direction. At fiqcD = 6.0, 
the finite temperature transition occurs somewhere in between temperature 
T = IjNt = 1/8 and r = 1/6 in lattice units. 

Each of the gauge configurations is independently fixed to Landau (and 
MA) gauge. As such, if the configurations are contaminated with Gribov 
copies [9] it is likely that the Gribov charge would vary between configura- 
tions [5]. Therefore, if Gribov copies have a substantial effect on Dj, their 
infiuence would show up in the size of our jackknife error bars depicted in 
Eigures 1-3. In fact, our errors are no bigger than errors for comparable 
gauge invariant quantities which are unaffected by Gribov copies. As it is, 
we observe no evidence for any anomalous Dj fiuctuation which might be 
due to the presence of Gribov copies in a subset of configurations. 

We define fractal dimension Dj of defects operationally as follows. Choose 
a cutoff value -Fmin and designate all sites x such that F^{x) < -Fmin as resis- 
tant. Let denote the sum over all resistant sites d and Nd{r) the number 
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of resistant sites inside a D = 3 + 1 dimensional hypersphere of radius r 
centered at resistant site d. Then the fraction of resistant sites inside the 
whole collection of such 3 + 1 dimensional hyperspheres of radius r is 



(4) 



d d 



The fractal dimension Dj of defects can be extracted as 




dlog N{r) 
J log r 



(5) 



Expression (5) obviously gives the right average value for the fractal di- 
mensionality if the fractals are sufficiently dilute and r is smaller than the 
size of the fractals. As r becomes much larger than the mean fractal separa- 



confining Nf = 40 lattices for both the Landau and the MA gauge defects. 
For comparison, the Dj of pseudo-defects comprised of randomly placed, ho- 
mogeneously distributed sites is also shown. Brackets ( ) indicate an average 
over importance sampling QCD gauge configurations. Only Landau gauge 
reveals a (nearly) fiat nonzero plateau region at small r, that is, only Lan- 
dau gauge exhibits definite clustering. In contrast MA gauge defects and, as 
required, the pseudo-defects do not exhibit a plateau. 

As illustrated in Figure 1, log(r/a) = 0.41 is well inside the r-plateau re- 
gion of Df for Landau gauge defects. Figure 2 depicts Dj at log(r/a) = 0.41 
as a function of /?, the number per unit volume of resistant sites, for Landau 
gauge defects at four temperatures and the homogeneous distribution. (For 
Landau gauge p is varied by adjusting cutoff -Fmin-) In contrast to the ho- 
mogeneous distribution where Dj monotonically rises with /?, Landau gauge 
defects all exhibit a plateau over a range of p. 



tion, Df approaches 4. Figure 1 shows 



dlogr 



function of log(r) on the 
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As indicated in Figure 2, the same temperature dependence is quali- 
tatively reproduced over a wide range of p. In Figure 3 we display Dj at 
log(r/a) = 0.41 and p ~ .005 for a range of temperatures in Landau gauge 
and MA gauge. As shown in Figure 3, only Dj in Landau gauge rises dra- 
matically as the temperature is lowered below the critical temperature. This 
drop correlates the clustering of resistant sites in Landau gauge with the 
onset of confinement. In contrast, MA gauge defects have Dj < 1 at all 
temperatures. 

We look for the aforementioned possible spacetime correlation between 
the ~ 1 Landau gauge defects and the Abelian monopole currents 
in MA gauge as follows. At each dual site x where k^[x) ^ we draw a 
3-dimensional sphere of radius r in the 3-surface orthogonal to p. Call this 
collection of spheres B{r). Counting the fraction of Landau gauge resistant 
sites inside B(r), we find that for p ~ .005 less than 10% of the resistant sites 
are contained inside B{r) even if r = 3 lattice spacings. Thus, Landau gauge 
defects are not spatially coincident with MA gauge monopole currents. As 
previously discussed, this does not necessarily mean they are not manifesta- 
tions of the same underlying QCD property in two different guises. At this 
point, we do not know. 

In summary, we have shown that the Landau gauge vacuum of lattice 
QCD has string-like defects. Like the Dirac string-originated defects of com- 
pact QFD, we suspect these QCD defects indicate the existence of an un- 
derlying gauge variant string structure in QCD, be it monopole currents or 
something else. It will be important to investigate what effects these defects 
have on QCD correlation functions. Based on the compact QFD results of 
Ref. [4], it is attractive to speculate that such strings may be (partially) 
responsible for disordering gluon and quark propagators giving rise to the 
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calculated nonzero gluon and quark mass poles in QCD [10]. 
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— pa =.005 24 x40 OLandau,F XMA,F ^l^homogeneous distribution — 




Figure 1: Landau and MA gauge are numerically achieved by maximizing 
T and T respectively with the same checkerboard relaxation procedure. 
The Landau gauge curve with resistant sites defined by shows & Dj ^ 
plateau over four r-slices before beginning to cross over to = 4. In 
contrast, both the MA gauge curve and the homogeneous distribution curve, 
whose resistant sites are randomly generated with a homogeneous probability, 
starts crossing over immediately from ~ at small r to 4 at large r with 
no plateau region, "a" refers to the lattice spacing at /3 = 6.0. Average ( ) 
is over Monte Carlo gauge configurations and all errors are jackknife errors. 
The guide-to-eye line indicates the plaieau region. 
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Figure 2: When p is too small, Dj ~ in all cases. As p is tuned up- 
wards, the Landau gauge Dj approaches a nonzero plateau which depends 
on temperature. These Landau gauge plateau are consistent with clustering. 
The corresponding MA gauge curves (not depicted) look like the homogenous 
curve at all temperatures. The jackknife errors on the higher temperature 
data is larger than the Nt = 40 data because the higher temperature lattices 
contain fewer sites and, hence, worst statistics. 
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Figure 3: In these fiqcD = 6.0 SU (3) lattices the finite temperature transition 
occurs somewhere between r = .125 and r = .16 — the region not traversed 
by the guide-to-eye fines. As depicted, -D/ in Landau gauge rises sfiarply as 
T decreases from tfie large r deconfined pfiase to tfie small r confined phase. 
In contrast, -D/ in MA gauge is always small. The bold horizontal line is the 
value of Df for a homogeneous, random distribution. 
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